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SECTION  I 


INTRODUCTION 


During  the  past  decade,  the  numerical  solution  of  chemically  reacting 
fluid  dynamic  flow  field  problems  has  received  considerable  emphasis  for  a 
wide  variety  of  applications.  The  feasibility  of  such  calculations  has  been 
made  possible  by  the  increase  in  computation  speed  achieved  by  succeeding 
generations  of  computers.  Despite  this,  the  solution  of  multidimensional 
chemically  reacting  flows  imposes  severe  computational  constraints,  even 
for  the  fastec*-  machines  currently  available.  For  rate  limited  chemistry, 
the  primary  difficulty  is  due  to  restrictive  numerical  stability  requirements 
dictated  by  the  reaction  rates.  Progress  in  this  area  is  determined  by  ad- 
vances in  numerical  techniques  and  is  not  the  province  of  this  report. 

For  equilibrium  chemistry,  stability  requirements  are  dictated  by  the 
fluid  dynamics,  and  the  computational  penalty  arises  from  the  determination 
of  the  composition  of  the  gas  mixture.  That  is,  solution  of  the  elemental 
conservation  equations  yields  the  element  mass  fractions.  To  determine  the 
thermodynamic  and  transport  properties  of  the  gas  mixture,  the  species 
mass  fractions  must  be  known.  The  determination  of  the  equilibrium  compo- 
sition by  species  from  given  element  mass  fractions  is  usually  carried  out  by 
minimizing  the  Gibb's  free  energy  which  involves  a time  consuming  iterative 
solution.  It  is  the  purpose  of  this  report  to  present  a computationally  fast, 
direct,  noniterative  method  for  determining  the  equilibrium  composition  of  a 
gas  mixture  for  a commonly  encountered  three-element  system  (carbon- 
nitrogen- oxygen). 


Direct,  noniterative  methods  have  been  developed  for  the  equilibrium 
composition  of  nitrogen-oxygen  systems  (i.e.,  air).  References  1,  2,  and  3 
utilize  various  levels  of  approximation  and  choices  of  independent  variables  to 
obtain  the  equilibrium  properties  of  air  with  excellent  accuracy.  Faster 
I iteration  schemes  have  been  developed  for  other  elemental  systems  (various 

i combinations  of  H,  C,  N and  O)  by  algebraically  combining  a set  of  equations 

into  a single  equation  (Refs.  4-6).  This  approach,  while  considerably  faster 
than  more  general  methods,  is  still  much  slower  than  a direct,  noniterative 
scheme.  Furthermore,  it  requires  substantial  algebraic  manipulation  to 
obtain  the  iterative  equation  and  an  accurate  initial  guess  to  insure  conver- 
gence. To  the  author's  knowledge,  there  do  not  exist  any  direct,  noniterative 
' methods  for  a carbon-nitrogen-oxygen  system.  This  report  provides  such  a 

I method. 

I 

I 
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States  Section,  The  Combustion  Institute,  1959. 
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The  solution  of  fluid  dynamic  problems  with  equilibrium  chemistry 
generally  falls  into  two  different  categories,  depending  on  the  interest  of  the 
investigator.  For  many  applications,  the  details  of  the  chemistry  are  of 
prime  importance,  particularly  with  regard  to  the  concentration  of  certain 
trace  species.  In  such  cases,  the  specification  of  the  equilibrium  composi- 
tion must  be  exact,  and  there  is  no  option  but  to  consider  a large  number  of 
species  and  solve  the  system  exactly  (using  Gibb's  free  energy  minimization 
technique,  for  example).  On  the  other  hand,  there  exist  many  problems  in 
which  the  main  interest  is  the  fluid  dynamic  quantities,  and  the  chemistry  is 
important  only  to  the  extent  that  it  determines  the  thermodynamic  and  trans- 
port properties  of  the  ga,s  mixture.  In  this  situation,  trace  species  may  be 
ignored  since  they  exert  no  significant  effect  on  the  thermodynamic  quantities. 
Furthermore,  moderate  errors  in  the  species  concentrations  may  be  tolerated 
since  the  mixture  properties  tend  to  average  the  errors.  This  is  particularly 
true  for  systems  in  which  the  elements  are  adjacent  to  one  another  in  the 
periodic  table  (such  as  carbon,  nitrogen,  and  oxygen).  This  report  presents 
an  approximate  chemistry  model  suitable  for  use  with  this  second  problem 
category. 

To  achieve  the  desired  objectives  (simplicity  and  fast  computation  time), 
two  assumptions  have  been  made,  consistent  with  the  intended  usage  of  the 
method.  That  is,  neither  ionized  nor  condensed  species  are  considered. 
Furthermore,  species  which  occur  in  trace  amounts  are  ignored.  The  major 
constraint  imposed  on  the  formulation  was  the  requirement  that  the  resulting 
equations  yield  analytic  solutions.  Since  analytic  solutions  do  not  exist  for 
polynomials  of  fifth  degree  and  higher,  the  problem  was  formulated  such  that 
only  cubics  and  quartics  are  encountered.  This  was  achieved  by  dividing  the 
problem  into  carbon-rich  and  oxygen-rich  domains  and  treating  the  dominant 
reactions  sequentially  rather  than  simultaneously  (i.  e.,  exact). 
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The  choice  of  independent  variables  is  dependent,  to  some  extent,  on 
the  intended  usage  of  the  method.  Since  the  partition  functions,  equilibrium 
constants  and  thermodynamic  properties  for  each  species  are  functions  of 
temperature,  the  temperature  must  be  one  of  the  independent  variables.  The 
other  independent  variable  mav  be  pressure  or  density.  For  fluid  dynamic 
problems  in  which  the  pressure  is  given  (i.e.,  boundary  layer  problems)  or 
treated  as  one  of  the  independent  variables,  it  is  the  obvious  choice  as  the 
variable  for  the  equilibrium  chemistry  calculation.  The  density  is  the  logical 
choice  of  variable  in  situations  where  the  continuity  equation  is  solved  for 
density  (i.e.,  flow  field  calculations).  | 

Since  the  choice  of  variables  exerts  a tremendous  influence  on  the 
complexity  of  the  model,  the  formulation  is  divided  into  two  parts  depending  ‘■ 

upon  whether  pressure  or  density  is  used  as  one  of  the  independent  variables. 

The  formulation  using  pressure  and  temperature  as  the  variables  is  the  most 
complex,  and  is  considered  first.  Following  the  development  of  the  model,  i 

extensive  comparisons  are  made  with  exact  solutions  to  assess  the  accuracy 
of  the  method. 
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SECTION  II 


PRESSURE  AND  TEMPERATURE  AS  INDEPENDENT  VARIABLES 


As  discussed  previously,  the  objective  of  this  report  is  to  calculate  the 
equilibrium  species  mass  fractions  and  mixture  properties  for  a carbon- 
nitrogen-oxygen  system.  For  general  use,  an  inert  monatomic  species  is 
included  in  the  formulation.  Since  the  inert  species  is  nonreactive  and  ioni- 
zation is  ignored,  it  is  characterized  solely  by  its  molecular  weight.  The 
pressure,  temperature  and  elemental  mass  fractions  are  assumed  to  be 
specified.  By  definition,  the  elemental  mass  fractions  (denoted  by  a super- 
script tilde)  must  sum  to  unity.  That  is 


* K„  + Kq  + iCj  = 1 


(1) 


where  the  subscript  I denotes  the  inert  species. 

The  critical  assumption  of  the  model  is  that  the  problem  may  be  divided 
into  two  distinct  regimes  (denoted  as  oxygen-rich  and  carbon- rich)  which  are 
essentially  uncoupled  from  one  another.  This  division  is  based  on  the  obser- 
vation that  carbon  will  always  react  preferentially  with  oxygen  to  form  CO 
before  reacting  with  nitrogen.  Therefore,  if  sufficient  oxygen  is  present  in 
the  system,  such  that  all  the  carbon  is  bound  up  in  CO  (i.  e.,  oxygen-rich 
system),  then  only  trace  amounts  of  carbon-nitrogen  compounds  (CN,  C2N2) 
and  carbon  sublimation  species  (C^,  C^)  will  be  present.  However,  if  suffi- 
cient carbon  is  present  in  the  system,  such  that  all  the  oxygen  is  bound  up  in 
CO  (i.e.,  carbon-rich  system),  the  excess  carbon  will  form  carbon  sublima- 
tion species  and  carbon-nitrogen  species.  Therefore,  a parameter  (denoted 
as  £)  defines  the  two  regimes 

^ “ 4 ^O  " ^C 


Oxygen- rich  system:  € > 0 
Carbon- rich  system:  € < 0 


(2) 


T 


As  mentioned,  the  two  regimes  are  uncoupled  from  one  another  and 
are  considered  separately  below. 

A.  OXYGEN-RICH  SYSTEM;  € >0 

In  view  of  the  above  comments,  the  carbon-nitrogen  and  carbon  subli- 
mation species  are  ignored.  Free  monatomic  carbon  does  not  exist  until  CO 
dissociates.  In  order  to  obtain  analytic  solutions,  the  dominant  reactions 
must  be  uncoupled  from  one  another.  The  CO2  and  O2  dissociation  reactions 
occur  at  low  to  moderate  temperatures  and  are  essentially  completed  before 
the  onset  of  N2  and  CO  dissociation  at  high  temperatures.  Therefore,  the 
reactions  will  be  treated  sequentially  in  pairs. 

1 . CO2  and  O2  Reactions 

The  reactions  are  denoted  as 


2 CO2— ►Z  CO  + O2 


O2-2O 


(3) 


While  these  reactions  are  taking  place,  all  nitrogen  is  presumed  to  be 
bound  up  in  N2.  The  elemental  mass  balances  yield 


12 

= 2F 

^o 

II 

0^ 

K,, 

N 

N. 

^1 

II 

12 


16 


32 


(4) 


The  molecular  weight  is  given  by 


1 , '"°z  , *^0  , '"co  , 

M " 28  32  IF"  28 


44 


(5) 
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The  pressure  equilibrium  constants  for  the  reactions  are  used  to 
relate  the  constituent  mass  fractions  of  the  products  and  reactants. 


^€©2  ""  ®1^  ^©2^00 


(6) 


a.  s 


121  P 


1 1568  K, 


Ot  s ^ 

2 8 


The  above  system  may  be  solved  to  obtain  the  four  species  mass  frac- 
tions and  molecular  weight  (M)  for  given  values  of  the  element  mass  fractions, 
pressure  and  temperature.  The  temperature  is  required  to  evaluate  the 
pressure  equilibrium  constants  (Kp^  and  Kp^).  The  pressure  equilibrium 
constants  are  presented  in  Appendix  A. 


The  solution  is  discussed  in  Appendix  B and  is  obtained  in  terms  of  a 
cubic.  Although  the  solution  yields  values  for  all  the  mass  fractions  of 
Eq.  (4),  only  the  values  of  CO2,  and  CO  (at  low  temperatures)  are  the 
final  values.  The  mass  fraction  of  atomic  oxygen  and  the  molecular  weight 
are  only  temporary  values  which  are  used  to  obtain  these  final  values.  The 
asymptotic  behavior  of  the  above  system,  for  both  large  and  small  tempera- 
tures, is  obtained  as 


For  T — ►0: 


Kq  = 0 


K 


- K iic 
O2  ■ O ■ 3 


(7) 


K 


11  ~ 
— — \c 

CO2  3 
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It  is  apparent  from  Eq.  (7)  that  the  mass  fraction  of  molecular  oxygen 
will  be  negative  for  sufficiently  low  temperatures  if  Kq  < 3Kq/8.  This  indi- 
cates that  insufficient  oxygen  exists  to  tie  up  all  the  carbon  in  In  reality, 

all  the  oxygen  would  combine  with  carbon  to  form  the  maximum  amount  of 
CO2,  and  the  remaining  carbon  would  condense  as  solid  carbon.  Thus,  even 
though  the  system  is  oxygen- rich  with  respect  to  CO  at  moderate  to  high 
temperatures,  it  is  also  carbon-rich  with  respect  to  CO2  at  low  temperatures. 
However,  since  this  occurs  at  temperatures  below  the  range  of  interest  of 
this  report,  it  causes  no  problem. 

As  the  temperature  approaches  infinity,  the  mass  fractions  for  O2  and 
CO^  have  the  correct  behavior.  The  carbon  is  all  contained  in  CO,  and  its 
mass  fraction  approaches  the  maximum  possible  value.  In  actuality,  at  high 
temperatures,  dissociation  would  take  place.  The  carbon  dioxide  dissocia- 
tion reaction  is  examined  in  Subsection  A-2. 

The  only  remaining  limit  of  interest  is  the  boundary  between  the  oxygen- 
rich  and  carbon- rich  domains  (i.e.,  € = 0).  The  above  system  has  the  follow- 
ing behavior  as  the  carbon-rich  domain  is  approached. 

£ -►0: 


0 

II 

0 

*^co 

0 

II 

0 

^co 

7 ~ 7 ~ 

= I^c*?^o 
= 0 


(9) 
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! i 
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i 


i 
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Again,  it  is  apparent  that  O2  and  CO^  are  consistent  with  the  postulated 
model,  namely,  that  all  the  oxygen  must  be  contained  in  CO  in  the  carbon- 
rich  regime.  Clearly,  at  sufficiently  low  temperatures,  this  will  be  incor- 
rect since  CO2  will  form  to  the  maximum  possible  extent  and  the  excess 
carbon  will  condense  as  solid  carbon.  However,  as  mentioned  earlier  this 
occurs  below  the  temperature  regime  of  interest. 

For  negative  values  of  €(i.e.,  carbon-rich),  the  system  yields  obviously 
incorrect  values  (negative  or  imaginary),  and  the  mass  fractions  of  O2  and 
CO2  are  set  equal  to  zero  without  solving  the  cubic. 

2.  N2  and  CO  Reactions 

The  reactions  considered  in  this  section  are  given  by 


N2-2N 


CO  ->c  + o 


(10) 


The  nitrogen  and  carbon  monoxide  dissociation  reactions  take  place  at 
high  enough  temperatures  such  that  molecular  oxygen  and  carbon  dioxide 
may  be  ignored.  Then  the  elemental  mass  balances,  molecular  weight  and 
equilibrium  relationships  are  given  by 


^0  - 

^0  28  *^CO 

- 

+ *^N2 

^C  = 

12 

K + •*  K 

C 28  ^CO 

II 

l-H 

1 

M = 

, '"”2  . KcO 

14  28  ' 28 

^CO 

= «3M 

= a.M  K^- 
4 N 

“3  • 

7P 

48  K 

^3 

^C 

12 


'TT 


M, 


(11) 


^4  = 7 K, 
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The  solution  to  this  system  is  obtained  in  terms  of  a quartic  and  is 
detailed  in  Appendix  C.  Since  the  resultant  equation  is  a quartic,  more  than 
one  solution  branch  is  possible.  This  is  indeed  the  case,  with  two  solution 
branches  for  the  oxygen- rich  and  carbon- rich  situations.  This  is  most  easily 
seen  from  the  asymptotic  behavior  of  the  quartic  as  the  temperature  vanishes 
(see  Appendix  C).  This  yields 

T-^0: 


€ >0  (oxygen- rich): 


= 0 


^CC)  ' 3 


Kc  = 0 


esi  4-  ^ 


i : 


’1 

1 


Since  the  quartic  possesses  valid  solution  branches  for  both  the  carbon- 
rich  and  oxygen-rich  case,  it  must  be  solved  for  either  case  to  obtain  the 
carbon  monoxide  mass  fraction. 

For  the  oxygen- rich  case,  a comparison  of  Eq.  (12)  with  Eq.  (8)  shows 
that  the  mass  fraction  of  CO  has  the  correct  asymptotic  behavior.  That  is, 
the  high  temperature  limit  of  Eq.  (8)  overlaps  the  low  temperature  limit  of 
Eq.  (12).  Thus,  the  final  value  for  the  CO  mass  fraction  (for  € > 0)  is  the 
smaller  of  the  two  values  obtained  from  Eqs.  (4)  and  (5)  and  Eq.  (11).  For 
the  carbon-rich  case,  the  final  value  for  CO  always  comes  from  11. 

The  behavior  of  the  system  at  large  temperatures  is  obtained  as 
T — ►oo: 

*^C  =^C 
*^0=^0 

Kn  = 0 
^CO  = 0 

At  high  temperatures,  there  is  only  one  valid  solution  branch  for  the 
system.  It  is  apparent  from  Eqs.  (12)  and  (13)  that  the  molecular  and  atomic 
nitrogen  mass  fractions  possess  the  correct  limiting  behavior  at  high  and 
low  temperatures . If  carbon-nitrogen  reactions  are  ignored,  these  are  the 
final  values  for  and  . If  carbon-nitrogen  reactions  are  included,  the 
final  values  for  the  nitrogen  mass  fraction  are  obtained  from  Subsection  B-1. 
This  is  discussed  further  in  Subsection  C. 


B.  CARBON-RICH  SYSTEM;  €<0 

It  was  seen  in  Subsection  A-1  that  the  and  CO^  vanished  identically 
as  the  boundary  between  the  carbon-rich  and  oxygen- rich  domains  (i.e., 

C = 0)  was  approached.  Thus,  for  negative  values  of  C,  the  mass  fractions  of 


-15- 
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©2  and  CO^  may  be  zeroed  and  the  cubic  need  not  be  solved.  This  yields 

values  for  O-  and  CO_  which  are  continuous  for  all  values  of  £. 

/i  c, 

j 

Since  the  carbon-nitrogen  and  carbon  sublimation  species  (as  calculated 
in  Subsections  B-1  through  B-3)  do  not  vanish  as  £ approaches  zero,  the 
converse  procedure  for  these  species  is  not  possible.  In  other  words,  zeroing 
these  species  for  positive  £ will  yield  discontinuous  values  at  £ = 0,  since  the  j 

solutions  do  not  yield  zero  values  as  the  boundary  is  approached  from  the 
carbon-rich  side.  Therefore,  in  order  to  obtain  continuous  solutions  for  all 
£,  the  carbon- rich  equations  of  this  section  will  be  solved  for  both  positive 
I and  negative  values  of  £.  That  is,  the  carbon- rich  equations  will  also  be 

solved  in  the  oxygen-rich  domain.  It  is  emphasized  that  the  values  obtained 
for  the  carbon-nitrogen  and  carbon  sublimation  species  (for  £ > 0)  must  be 
quite  small  since,  in  reality,  they  exist  in  only  trace  amounts  in  the  oxygen- 
rich  domain.  This  will  be  assured  by  utilizing  the  mass  fraction  of  CO2  (a 
known  quantity  from  Subsection  A-1)  in  the  elemental  carbon  mass  balance. 

As  mentioned  in  Subsection  A-2,  the  set  of  Eqs.  (11)  is  also  solved  for 
the  carbon- rich  case  to  obtain  the  carbon  monoxide  and  atomic  oxygen  mass 
fractions.  After  obtaining  these  values,  they  are  then  treated  as  known 
quantities  in  the  carbon  and  oxygen  elemental  mass  balances  for  the  remainder 
of  the  carbon- rich  calculations.  If  the  carbon-nitrogen  reactions  are  ignored, 
then  the  values  for  atomic  and  molecular  nitrogen  mass  fractions  are  obtained 
from  Subsection  A-2.  Generally,  however,  the  dominant  carbon-nitrogen 
reactions  are  included  and  are  used  to  obtain  the  nitrogen  mass  fraction. 

It  was  noted  previously  that  the  formation  of  solid  carbon  was  possible 
for  the  oxygen-rich  system  at  low  temperatures  where  the  carbon  was  bound 
up  in  CO2  instead  of  CO,  For  the  carbon-rich  situation,  the  formation  of 
solid  carbon  is  possible  at  much  higher  temperatures.  At  sufficiently  high 
temperatures,  the  excess  carbon  is  contained  in  the  carbon  sublimation 
species  (C2  and  C^).  However,  as  the  temperature  is  lowered  these  species 
vanish  and  the  excess  carbon  condenses  out  as  solid  carbon.  The  tempera- 
ture at  which  this  occurs  is  dependent  upon  the  pressure,  and  sets  a lower 
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limit  for  the  validity  of  the  present  model.  This  lower  linx...  .cmperature  is 
always  high  enough  so  that  molecular  oxygen  and  carbon  dioxide  exist  only  in 
trace  amounts.  This  is  consistent  with  zeroing  these  mass  fractions  in  the 
carbon-rich  domain. 

In  order  to  obtain  analytic  solutions,  the  and  reactions  must  be 
uncoupled  from  the  carbon-nitrogen  reactions.  The  sequence  in  which  the 
reactions  are  considered  follows  from  the  observation  that  at  low  carbon 
mass  fractions  CN  is  predominant,  whereas  and  are  predominant  for 
large  carbon  mass  fractions.  It  is  presumed,  therefore,  that  the  carbon- 
nitrogen  reactions  have  first  call  on  the  available  carbon  and  will  form  CN 
to  the  maximum  possible  amount,  subject  to  equilibrium  constraints.  Follow- 
ing the  carbon-nitrogen  reactions,  the  remaining  available  carbon  will  be 
absorbed  in  C2  and  C^. 

1 . N;>  and  CN  Reactions 

The  dominant  carbon-nitrogen  reaction  is  represented  by  the  following 
reactions: 

N^^ZN 

C^— ►2C  (14) 

CN-*C  + N 

The  elemental  mass  balances,  molecular  weight  and  equilibrium  rela- 
tionships are  given  by 

^N  “ ^ *^N2  ZS"  ^CN 


K.,  ^N,  K K 

^ 4. ?.  j. + — — + — — + — L -I-  B 

14  ^ 28  26  12  24  ^ Mj  ^2 


^^2 

= a.M 

4 N 

= a.M 

^2 

5 C 

= Ot/M  K^K 

CN 

6 C 

“5  " 

P 

6Kp 

^*6  “ 84  K, 


As  indicated,  the  functions  and  are  obtained  from  the  solutions  of 
Subsections  A-1  and  A-2  and  are  given  by 


8 = 11k  +11k 

^1  “ 28  CO  44  '^CO, 


^CO  ^‘=°2 

^2  = 28  TT'  32  44 


The  solution  to  this  system  is  in  terms  of  a quartic  and  is  contained  in 
Appendix  D.  The  limiting  behavior  for  large  and  small  temperatures  is 
given  by 

T-^0: 


Kc  = 0 

Kn  = 0 
^CN  = 0 


- ^1 

= I Kq  for  € < 0 
ySj  = for  £ > 0 


(17  continued) 
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T-^c»: 


^N  = 
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0 

II 

^CN  = 
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II 

u 
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j 

II 

0 

It  is  apparent  that  at  low  temperatures,  for  a carbon-rich  system,  the 
excess  carbon  is  absorbed  in  instead  of  condensing  as  solid  carbon.  If 
the  system  is  oxygen-rich,  should  vanish  with  the  temperature  and  it  does. 
The  results  for  nitrogen  and  cyano  (CN)  have  the  correct  limiting  behavior  at 
both  temperature  extremes.  The  final  results  for  the  mass  fractions  of 
molecular  and  atomic  nitrogen  are  obtained  from  this  system.  The  results 
for  CN,  and  C are  temporary  values  which  are  used  only  to  obtain  the 
result.  The  value  obtained  for  the  CN  mass  fraction  is  quite  accurate,  except 
in  situations  where  cyanogen  (C2N2)  forms  in  significant  amounts.  Therefore, 
in  the  following  sections  the  nitrogen  bound  up  in  CN  is  reapportioned  between 
CN  and  C2N2. 

2.  C2N2  and  CN  Reactions 

Cyanogen  will  form  in  moderate  amounts  close  to  the  temperature 
boundary  where  condensed  carbon  will  form.  In  this  subsection,  a small 
correction  is  applied  to  the  mass  fraction  of  cyano  obtained  in  Subsection  B-1. 
The  mass  fractions  of  atomic  and  molecular  nitrogen  are  known  from  Sub- 
section B-1.  In  the  interest  of  simplicity,  the  molecular  weight  of  Subsection 
B-1  is  also  used  for  this  calculation.  The  reaction  of  interest  is  given  by 

C^N^-^Z  CN  (18) 


The  nitrogen  mass  balance  and  the  equilibrivim  relationship  are  given  by 


” 25"  ^CN  52  ^C2N2  ^ 
^^2^2  ^ 


P 

^7  “ 13  K, 


^3  ■'■  *^N- 


The  system  yields  the  following  quadratic 


^13  “7^^  *^CN  13  ^CN  " ■ ^3  j ° 


Since  M is  taken  from  the  result  of  Subsection  B-1,  this  equation  may  be 
solved  for  the  CN  mass  fraction.  The  cyanogen  mass  fraction  is  obtained 
from  the  equilibrium  relationship.  These  are  the  final  values  for  CN  and 
C2N2.  Since  approaches  Kj^  at  both  temperature  extremes,  it  is  easily 
seen  that  the  cyano  and  cyanogen  mass  fractions  vanish  for  large  and  small 
temperatures. 

3.  C2  and  C-^  Reactions 

In  this  subsection,  the  remaining  available  carbon  is  apportioned 
between  C,  C2  and  C^.  The  pertinent  reactions  are 

C--»2C 

^ (21) 

C3-»3C 


The  system  is  determined  from  the  following  set  of  equations: 
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The  P functions  are  known  from  the  preceding  calculations  and  are 
given  by 
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If  the  carbon-nitrogen  reactions  are  ignored,  the  mass  fractions  for 
CN  and  C2N2  are  zeroed  in  the  above  expressions  for  the  fi  functions. 

The  solution  is  obtained  in  terms  of  a cubic  and  is  presented  in 
Appendix  E.  The  limiting  behavior  for  the  system  is  given  by 
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K_  =0 

^€3  = - ^4 

^4  = IV 
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Thus,  at  low  temperatures,  for  a carbon-rich  system,  the  excess  car- 
bon is  contained  in  C3  instead  of  condensing  out  in  the  form  of  solid  carbon. 
For  an  oxygen-rich  system,  the  C3  vanishes  with  the  temperature,  consistent 
with  the  model. 

If  the  quantity  is  negative,  it  indicates  that  all  the  carbon  is 

bound  up  in  carbon  monoxide,  carbon  dioxide  and  the  carbon  nitrogen  com- 
pounds, and  none  is  available  for  the  formation  of  and  C3.  In  this  case, 
the  mass  fractions  for  C2  and  C3  are  zeroed,  and  the  cubic  is  not  solved. 

C.  COMPUTATIONAL  SEQUENCE 

It  is  important  to  note  that  the  computational  sequence  of  Subsections  A 
and  B is  not  unique  and  may  be  changed  to  reflect  different  priorities.  For 
example,  the  calculation  order  of  these  subsections  was  chosen  to  insure 
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continuity  at  the  dividing  line  between  the  oxygen- rich  and  carbon- rich 
domains.  Before  discussing  this  further,  the  computation  sequence  is  sum- 
marized as  follows. 

For  positive  C,  the  mass  fractions  of  and  CO^  are  obtained  from 
Subsection  A-1.  For  negative  £,  these  values  are  zeroed.  For  positive  €, 
the  CO  mass  fraction  is  the  smaller  of  the  two  values  from  Subsections  A-1 
and  A-2.  For  negative  €,  the  CO  value  comes  from  Subsection  A-2.  In  the 
temperature  regime  of  interest  (high  enough  so  that  condensed  carbon  is  not 
present)  the  CO  mass  fraction  will  be  continuous  across  the  boundary  C a 0. 

If  the  carbon-nitrogen  reactions  are  ignored,  the  mass  fractions  for  N and 
come  from  Subsection  A-2  for  all  €.  If  the  reactions  are  included,  the 
nitrogen  values  come  from  Subsection  B-1  for  all  €.  The  values  for  CN  and 
C^N^  are  solved  for  both  positive  and  negative  € from  Subsection  B-2.  The 
mass  fractions  of  and  C^  are  also  solved  for  all  € and  are  obtained  from 
Subsection  B-3.  The  atomic  oxygen  and  carbon  mass  fractions  are  obtained 
from  Eqs.  (25). 

As  mentioned  earlier,  the  equations  for  the  carbon- rich  domain  are 
also  solved  in  the  oxygen- rich  region  to  insure  continuity  across  the  boundary 
£ = 0.  If,  however,  the  problem  of  interest  always  lies  in  the  oxygen- rich 
domain,  such  that  the  dividing  line  between  oxygen  and  carbon  rich  is  never 
crossed,  then  the  carbon- rich  calculations  would  be  skipped.  The  nitrogen 
mass  fractions  would  then  be  obtained  from  Subsection  A-2.  For  example, 
any  fluid  dynamics  problem  with  a fixed  set  of  elemental  mass  fractions  (such 
that  £ > 0)  would  fall  into  this  category. 

D.  MIXTURE  PROPERTIES 

All  the  molecular  mass  fractions  were  computed  in  Subsections  A and  B. 
Following  these  calculations,  the  complete  elemental  mass  balances  are  used 
to  compute  the  atomic  mass  fractions. 
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The  final  value  for  the  mixture  molecular  weight  is  computed  from 

K. 


M ~Z^  M. 


(26) 


The  mixture  gas  constant  is  obtained  from 

R 


R = 


M 


(27) 


where  R^  is  the  universal  gas  constant. 

The  mixture  enthalpy  is  calculated  from 
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The  enthalpy  of  each  species  is  obtained  from  the  partition  functions 
of  Appendix  A,  That  is 
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The  equation  of  state  relates  the  pressure  and  density. 

PM  = R^PT  (30) 

Some  discussion  as  to  usage  is  in  order.  As  indicated  earlier,  one  of 
the  independent  variables  must  be  temperature  because  the  thermodynamic 
data  for  the  individual  species  are  specified  in  terms  of  temperature.  If  the 
energy  equation  is  solved  for  temperature,  then  it  is  immediately  available 
for  input  into  the  chemistry  calculation.  However,  the  energy  equation  is 
usually  solved  for  enthalpy  and,  thus,  the  temperature  is  not  known.  To 
obtain  the  temperature  for  input  into  the  chemistry  subroutine,  the  enthalpy 
at  a particular  mesh  point  is  expanded  about  the  last  iterate  level,  and  the 
resulting  Taylor  series  is  solved  for  temperature. 

That  is,  the  enthalpy  is  specified  as 

H = H{T,P,K.).  (31) 

At  any  mesh  point  in  the  flow  field,  this  expression  may  be  expanded  in  a 
Taylor  series  about  the  last  iterate  (or  time)  level.  If  the  pressure  is 
assumed  to  be  fixed  and  does  not  vary  with  the  iterate  level  (i.e.,  boundary 
layer  problem),  the  expansion  is  obtained  as 
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The  iteration  levels  are  denoted  by  n for  the  last  iteration  and  n+1  for 
the  current  iteration.  The  enthalpy  and  elemental  mass  fractions  are  available 
at  the  current  iteration  from  the  energy  and  elemental  conservation  equations, 
and  all  quantities  are  known  at  the  last  iteration.  Therefore,  the  equation 
may  be  solved  for  temperature  at  the  current  iteration  and  input  to  the  chem- 
istry calculation. 

The  indicated  derivatives  (Cp  and  Z)  are  evaluated  as  follows: 


J'T,P 


(33) 


The  derivatives  may  be  obtained  by  differentiating  the  analytic  solution. 
It  is  important  to  note  that  these  derivatives  may  be  calculated  without  resolv- 
ing the  cubics  and  quartics.  Therefore,  evaluating  these  quantities  analyti- 
cally, rather  than  numerically,  will  substantially  decrease  the  computation 
time. 

E.  CLOSURE 

In  order  to  obtain  equations  which  were  of  fourth  order  or  lower,  it 
was  necessary  to  uncouple  the  reactions  and  consider  them  sequentially.  The 
motivation  for  determining  which  reactions  to  couple  and  treat  independently 
of  the  other  reactions  was  discussed  above.  A considerable  effort  was  made 
to  examine  alternate  possibilities  for  increased  accuracy  and  simplicity, 
several  of  which  are  described  here. 

An  attempt  was  made  to  consider  a single  set  of  equations  for  the 
oxygen- rich  regime  [consisting  of  Eqs.  (3)  and  (10)]  and  a single  set  for  the 
carbon-rich  regime  [Eqs.  (10),  (14),  (18)  and  (21)].  In  both  cases,  poly- 
rlomials  of  higher  than  fourth  order  were  encountered.  The  same  difficulty 
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was  encountered  when  trying  to  treat  the  carbon-nitrogen  reactions  [Eqs.  (14) 
and  (18)]  simultaneously.  An  attempt  to  incorporate  the  CN  reaction  of  i 

Eq.  (14)  with  the  CO  and  N2  reaction  of  Eq.  (10)  yielded  a fifth  order  poly- 
nomial. Therefore,  it  is  felt  that  the  present  model  is  as  simplified  as 
possible  within  the  fourth  order  polynomial  constraint  imposed  on  the  system. 

In  several  places  during  the  formulation,  allusions  were  made  to  the 
possibility  of  ignoring  the  carbon-nitrogen  reaction.  These  reactions  were 
included  for  generality  and  in  the  anticipation  that  they  might  prove  signifi- 
cant for  a carbon-nitrogen  system.  It  was  found,  however,  that  these  reac- 
tions exert  a negligible  effect  upon  the  mixture  properties  when  the  system  is  I 

nearly  all  carbon  and  nitrogen  with  small  amounts  of  oxygen.  For  a carbon- 
air  system,  with  a 50  percent  carbon  mass  fraction,  it  is  shown  in  Section  IV 
that  the  inclusion  of  the  carbon-nitrogen  reactions  has  only  a small  influence 
on  the  mixture  properties  (maximum  difference  of  two  to  three  percent).  An  j 

additional  set  of  comparisons  was  made  for  a carbon-nitrogen  system  with  ; 

the  same  conclusion.  It  is  believed,  therefore,  that  the  carbon-nitrogen 
reactions  may  be  ignored  with  complete  generality. 

Thus,  in  its  simplest  form  (ignoring  the  carbon-nitrogen  reaction),  the 
mixture  properties  of  a carbon-nitrogen-oxygen  system  may  be  determined 
by  solving  two  cubics  (Subsections  A-1  and  B-3)  and  a quartic  (Subsection  A-2). 

It  is  noted  that  the  formation  of  nitric  oxide  (NO)  was  ignored.  Hansen's  , 


results  (Ref.  1)  show  that  it  exerts  a negligible  influence  on  the  mixture  pro- 
perties of  air.  Thus,  since  the  nitrogen  must  compete  with  carbon  for  the 
available  oxygen  in  the  carbon-nitrogen-oxygen  system,  it  will  be  even  less 
significant  for  this  system.  However,  if  the  mass  fraction  of  NO  were  of 
interest,  it  could  be  determined  by  using  the  molecular  weight  and  CO  mass 
fraction  as  calculated  by  the  present  formulation  and  solving  the  elemental 
mass  balances  for  nitrogen  and  oxygen  (using  K^q  as  a known  quantity)  along 
with  the  equilibrium  relationships  for  the  N2,  O2  and  NO  dissociation  reac- 
tions . 
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In  the  Appendices,  an  implicit  assumption  is  made  that  nitrogen  was 
always  present  in  the  system  (i.  e.,  ^ 0).  Thus,  while  the  equations  are 

well  behaved  for  zero  values  of  the  carbon  or  oxygen  elemental  mass  frac- 
tions, the  equations  of  Appendix  C are  singular  for  a zero  nitrogen  mass 
fraction.  Furthermore,  if  = 0,  thj  CO  dissociation  reaction  is  treated  by 
itself  and  yields  a quadratic  equation  (see  Appendix  C). 

If  the  carbon  elemental  mass  fraction  is  zero,  the  present  results 
reduce  to  Hansen's  (Ref,  3)  model  (uncoupled  and  dissociation  reac- 
tions with  no  nitric  oxide).  Because  of  the  importance  of  the  nitrogen-oxygen 
I system  (which  contains  air  as  a special  case),  a more  complete  model  is 

treated  in  Appendix  F.  It  is  shown  that  the  N2,  and  NO  reactions  can  be 
treated  simultaneously  and  the  solution  obtained  in  terms  of  a quartic. 

Some  comments  concerning  round-off  errors  in  the  extraction  of  the 
roots  is  in  order.  Since  the  coefficients  vary  by  several  hundred  orders  of 
magnitude  over  the  temperature  range  of  interest  (500- 10,000°K),  the  equa- 
tions change  character  at  the  temperature  extremes.  The  quartics,  for 
example,  reduce  to  quadratics  for  large  and  small  temperatures.  Therefore, 
extracting  the  roots  from  the  analytic  solution  for  a quartic  introduces  some 
round-off  errors  in  this  situation.  After  obtaining  the  desired  root  from  the 
quartic,  it  was  tested  to  ensure  that  it  did  indeed  satisfy  the  equation  to 
within  some  tight  tolerance.  At  high  and  low  temperatures,  it  exceeded  the 
tolerance,  and  the  root  was  recalculated  using  the  asymptotic  forms  of  the 
equations  derived  in  the  Appendices. 

Since  the  analytic  solution  for  a quartic  is  only  marginally  more  com- 
plicated than  that  for  a cubic,  the  computation  time  to  obtain  the  four  roots 
is  negligible.  The  decision  as  to  which  of  the  real  roots  (imaginary  roots 
were  ignored)  is  the  desired  root  is  trivial.  In  all  cases,  the  smallest 
positive  root  is  the  correct  value. 
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SECTION  III  I 

DENSITY  AND  TEMPERATURE  AS  INDEPENDENT  VARIABLES  i 

i 

The  analysis  of  Section  II  is  repeated  here  with  the  density  replacing  j 

the  pressure  as  one  of  the  given  variables.  An  exannination  of  any  of  the  \ 

equilibrium  relationships  shows  that  the  product  of  pressure  and  molecular  | 

weight  always  appears  in  the  equation.  Thus,  for  a given  value  of  pressure,  j 

the  molecular  weight  remains  as  an  unknown  in  the  eqviilibrium  expression.  ,i 

This  couples  the  elemental  mass  balances  and  equilibrium  equations  to  the  | 

j 

equation  relating  molecular  weight  and  species  mass  fractions.  It  is  this  | 

coupling  of  equations  which  greatly  complicates  the  algebra  required  to  solve  | 

the  system.  i 

i 

When  the  density  and  temperature  are  given,  the  product  of  pressure  | 

and  molecular  weight  is  easily  evaluated  from  the  equation  of  state  [Eq.  (30)].  \ 

Therefore,  the  equilibrium  relationships  do  not  contain  the  molecular  weight  ' 

as  an  unknown  quantity,  and  the  system  may  be  solved  for  the  species  mass  ! 

fractions  independently  of  the  molecular  weight  eqviation.  Since  the  product  ' 

of  pressure  and  molecular  weight  is  a known  quantity,  it  is  absorbed  into  a I 

new  set  of  coefficients  denoted  as  3..  These  coefficients  are  simply  the  fit.  of  > 

1 1 j 

Section  II  multiplied  by  the  molecular  weight;  that  is  | 

S.  = Mof.  (34) 

The  partitioning  of  the  problem  into  carbon-rich  and  oxygen-rich 
regimes  and  the  order  in  which  the  reactions  are  considered  remains  un- 
changed. As  before,  the  oxygen-rich  reactions  are  considered  first. 

A.  OXYGEN- RICH  SYSTEM :€  > 0 

Due  to  the  uncoupling  of  the  elemental  mass  balance  equations  and 
equilibrium  expressions  from  the  molecular  weight  equation,  it  is  possible  to  ! 
treat  all  the  oxygen-rich  reactions  in  one  set  of  equations.  That  is,  the  CO 


dissociation  reaction  may  be  included  with  the  CO^  and  reactions  and  the 
solution  obtained  in  the  form  of  a quartic.  However,  the  negligible  increase 
in  accuracy  does  not  warrant  replacing  the  cubic  with  a quartic.  Therefore, 
the  CO  dissociation  reaction  is  treated  separately,  and  the  analysis  exactly 
parallels  that  of  Subsection  II-A. 

1.  CO2  and  O2  Reactions 

The  elemental  mass  balances  are  exactly  the  same  as  Eqs.  (4)  and  are 
not  repeated  here.  The  equilibrium  relationships,  with  the  newly  defined Ot., 
are  given  as 


= “2^ 


~ _ 121  PM 

1 " 1568  Kp 


~ _ PM 
«2=  8K 


Substitution  into  the  elemental  mass  balances  yields 


AKq  + + CKq  - I f = 0 


A = — 3 \/~3~a 
^ 112  12 


® ~ 7 “2  1 1 


_ _ 3 , / ' ~ ~ / 8 ~ 3 r; 

^ “ 7 “12  \ 1 1 ■ 1 1 


^CO  3 J_„ 

I^TT^O^  V2 


The  limiting  behavior  of  this  system  is  exactly  the  same  as  it  was  for 
pressure  as  the  independent  variable. 


2. 


N->  and  CO  Reactions 


When  the  pressure  was  used  as  the  variable,  the  N2  and  CO  reactions 
were  linked  together  only  through  the  molecular  weight  equation.  Since  the 
use  of  density  as  a variable  eliminates  the  molecular  weight,  these  two 
reactions  are  now  uncoupled  and  may  be  solved  separately. 

The  equilibrium  relationships  are  given  by 


^CO  = "3^0  ^O 


= 

4 N 


(37) 


where  the  Ct.  are  defined  according  to  Eq.  (34). 

Substitution  into  the  elemental  mass  balances  of  Eq.  (11)  yields  a set  of 
quadratic  equations. 


-f“3")  -^o  = ° 


(38) 


There  is  only  one  positive  root  for  Kj^,  and  this  is  the  correct  value. 
At  low  temperatures,  the  two  roots  for  K correspond  to  the  two  solution 
branches  for  the  oxygen-rich  case  (€  > 0)  and  the  carbon- rich  case  (€  < 0). 
The  oxygen-rich  case  corresponds  to  the  plus  sign  in  front  of  the  radical  in 
the  equation  for  the  roots  of  a quadratic.  The  carbon-rich  situation  cor- 
responds to  the  minus  sign.  At  high  temperatures  there  is  only  one  positive 
root  for  K^,  and  it  is  the  correct  value.  The  quadratic  for  Kq  has  similar 
behavior. 

The  asymptotic  behavior  of  this  system,  for  large  and  small  tempera- 
tures, is  identical  to  Eqs.  (12)  and  (13). 


B.  CARBON-RICH  SYSTEM;  e<  0 

Again,  the  carbon-rich  analysis  exactly  parallels  that  of  Subsection  II-B 
with  the  same  elemental  mass  balances  and  equilibrium  relationships  (with 
the  molecular  weight  absorbed  in  the  3^). 

1.  N2  and  CN  Reactions 

The  equilibrium  relationships  are 

K 


N2 


^C2  = 


(3  9) 


^CN  = "6  Vn 


Substitution  into  the  elemental  mass  balances  of  Eq.  (15)  yields 


AK^,  + BK^,  + CK^,  + DK-,  +3_K^-  = 0 
N N N N 5 N 


A « ^2  4Z  ^ 

^ = «5“4  - 1^«4«6 

23435-^3^-^348^ 

^ + S5  - 2S435Kj^  + 

■ ■n“6  " ll|“6  ^^C  • ^1^ 


(40) 


13 


The  carbon  mass  fraction  is  calculated  from  the  following  quadratic: 

(k^-^,)=0  (41) 


^54 


The  limiting  behavior  for  large  and  small  temperatures  is  given  by 
Eq.  (17). 


« 

fc 
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2.  C2N2  and  CN  Reactions 

The  purpose  of  this  calculation  is  simply  to  reapportion  the  nitrogen 
absorbed  in  CN  (in  Subsection  B-1)  between  CN  and  C^N^. 

Utilizing  the  nitrogen  mass  balance  and  the  equilibrium  relationship 

yields 


(w37)Kcn+T 


*^CN  13  *^CN  ■ ^^N  “ ^3^  ® ° 


The  single  positive  root  is  the  correct  value.  At  the  temperature 
extremes  both  the  CN  and  C2N2  mass  fractions  vanish. 

3.  C2  and  C3  Reactions 

The  equilibrium  expressions  are  given  by 

(43) 

^C3  = 4 

Substitution  into  the  carbon  mass  balance  of  Eq.  (22)  yields  the  fol- 
lowing cubic: 

- 0 <«) 

As  before,  if  ^ K^,  all  the  carbon  is  tied  up  in  CO  and  the  carbon- 
nitrogen  compounds  and  is  not  available  to  form  C2  or  low  tempera- 

ture, the  excess  carbon  is  contained  in  C^. 

C.  CLOSURE 

After  the  calculations  for  the  molecular  mass  fractions  are  completed, 
the  relationships  of  Subsection  II-D  are  used  to  calculate  the  atomic  mass 
fractions  and  mixture  properties. 
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The  computational  sequence  is  exactly  the  same  as  for  Section  II  and  is 
again  chosen  to  insure  continuity  across  the  boundary  € = 0.  The  calculation 
sequence  is  summarized  in  Subsection  II-C,  and  the  same  comments  con- 
cerning alternate  computational  sequences  are  applicable  here. 

Since  the  physical  model  treated  in  Sections  II  and  III  is  identical,  the 
behavior  of  the  system  was  not  discussed  in  Section  III.  It  should  be  empha- 
sized, however,  that  the  results  of  this  section  are  subject  to  the  same 
limitations  and  qualifications  as  those  of  Section  II.  In  particular,  the  same 
comments  apply  here  concerning  the  formulation  of  solid  carbon  and  the 
advisability  of  considering  carbon-nitrogen  reactions. 

A comparison  of  the  analysis  of  this  section  with  that  of  Section  II  (as 
contained  in  the  relevant  appendices)  makes  it  obvious  that  the  use  of  density, 
instead  of  pressure,  as  a variable  results  in  a tremendous  reduction  in  alge- 
braic complexity.  The  greatest  simplification  is  contained  in  Subsection  A-2 
where  the  quartic  of  Section  II  was  replaced  by  two  uncoupled  quadratic  equa- 
tions. Thus,  in  its  simplest  form  (ignoring  carbon-nitrogen  reactions)  the 
formulation  using  density  does  not  require  the  solution  of  any  quartics. 

It  is  noted  that  the  elimination  of  the  molecular  weight  from  the  deter- 
mination of  the  species  mass  fractions  means  that  the  inert  species  also  has 
no  influence  on  the  calculation  of  the  mass  fractions.  The  inert  species  mass 
fraction  is  required  only  to  calculate  the  molecular  weight  following  the  cal- 
culation of  the  other  mass  fractions. 


SECTION  IV 


RESULTS  AND  DISCUSSION 


In  this  section,  the  analytic  solution  developed  in  Section  II  (given 
pressure  and  tenaperature)  is  compared  with  exact  calculations  to  determine 
accuracy  and  computation  speed.  The  results  of  Section  III  (given  density 
and  pressure)  will  yield  exactly  the  same  accuracy  as  those  of  Section  II 
because  the  physical  model  is  identical  in  both  cases.  The  computation  speed 
for  given  density  will  be  somewhat  faster  than  for  given  pressure  because  of 
the  reduced  algebraic  complexity.  The  accuracy  of  the  analytic  solution  is 
considered  first. 

A.  ACCURACY 


There  exist  two  distinct  kinds  of  fluid  dynamic  problems  involving 
equilibrium  chemistry.  The  first  consists  of  the  situation  in  which  the  ele- 
mental mass  fractions  are  fixed  by  external  circumstances  and  remain 
constant  with  time  and  throughout  the  flow  field.  For  this  class  of  problem, 
a Mollier  diagram  may  be  generated,  for  a fixed  set  of  elemental  mass 
fractions,  to  provide  the  properties  of  the  gas  mixture  over  a range  of  tem- 
perature and  pressure.  Since  the  elemental  mass  fractions  are  known  and 
are  independent  of  the  fluid  dynamics  calculation,  the  Mollier  diagram  may 
be  precomputed  from  an  exact  equilibrium  chemistry  code  and  stored  on 
cards  or  magnetic  tape  for  use  with  the  fluid  dynamic  computation.  In  this 
manner,  the  chemistry  calculation  reduces  to  a two-dimensional  interpola- 
tion of  the  stored  data.  A commonly  encountered  problem  of  this  class 
consists  of  calculating  the  flow  field  about  a non-ablating  body  in  a planetary 
atmosphere  (or  in  a wind  tunnel,  ballistic  range  or  rocket  exhaust). 

In  the  second  problem  category,  the  elemental  mass  fractions  vary 
throughout  the  flow  field  and  are  determined  as  part  of  the  solution.  A com- 
mon example  of  this  involves  the  calculation  of  the  boundary  layer  over  an 
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ablating  surface  where  the  ablation  process  changes  the  elemental  mass 
fractions.  Since  the  amount  of  ablation  depends  on  the  heat  transfer  to  the 
surface,  it  is  apparent  that  the  determination  of  the  fluid  dynamic  and  ele- 
mental mass  fractions  is  a completely  coupled  process.  It  is  obviously  not 
possible  to  provide  computer  storage  for  separate  Mollier  diagrams  for  each 
possible  set  of  elemental  mass  fractions  that  might  be  encountered  on  a prob- 
lem of  this  kind.  In  this  case,  the  equilibrium  chemistry  calculations  must 
be  carried  out  simultaneously  with  the  fluid  dynamic  computation. 

Thus,  while  the  analytic  solution  of  this  report  may  be  used  with  either 
category  of  problem,  it  apparently  will  be  most  useful  for  those  problems 
with  varying  mass  fractions.  In  fact,  the  author's  prime  motivation  in 
developing  this  solution  was  for  its  use  in  boundary  layer  calculations  on  a 
carbon  surface  ablating  in  air. 

The  first  test  case  consists  of  a Mollier  diagram  for  a planetary  atmos- 
phere and  is  typical  of  the  first  kind  of  problem  (fixed  elemental  composition).  1 

This  case  was  denoted  as  a postulated  Mars  atmosphere  in  Ref.  7 and  con-  ^ 

sists  of  25  percent  N^,  32  percent  Ar  and  43  percent  CO^  by  volume.  The 
corresponding  elemental  mass  fractions  are 

= 0.  1333 
= 0.  1808 
Kq  = 0.3553 
Kj  = 0. 3306 

This  case  is  oxygen-rich  with  respect  to  the  formation  of  both  carbon  ' 

monoxide  and  carbon  dioxide.  That  is,  at  low  temperatures  all  the  carbon  is 
bound  up  in  CO2  and,  thus,  solid  carbon  will  never  be  present  in  this  mixture. 

The  results  for  molecular  weight  and  enthalpy  are  presented  in  Figs. 

1 and  2 for  a range  of  temperature  and  pressure.  The  solid  lines  are  the 
approximate  analytic  solution,  and  the  symbols  are  the  exact  result  from 

) 
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Ref.  7.  It  is  apparent  that  the  agreement  is  excellent,  except  at  low  pressure 
and  high  temperature  where  ionization  becomes  significant.  The  mass  frac- 
tions at  a pressure  of  ten  atmospheres  are  compared  in  Figs.  3a  and  3b.  The 
exact  calculations  for  the  mass  fractions  were  not  presented  in  Ref.  7 and 
were  obtained  from  the  NEST  (Ref.  8)  equilibrium  chemistry  code.  Again, 
the  agreement  is  excellent.  The  slight  error  in  the  mass  fractions  of  molecu- 
lar oxygen  and  nitrogen,  in  the  temperature  range  of  3000-6000°K,  is  due  to 
the  neglect  of  nitric  oxide  in  the  approximate  solution. 


The  remainder  of  the  cases  chosen  for  comparison  are  pertinent  to  the 
second  category  problem  (varying  elemental  mass  fractions)  discussed  above. 
The  cases  are  all  carbon-air  mixtures  with  varying  amounts  of  carbon.  The 
air  was  argon  free  and  the  volumetric  composition  was  80  percent  molecular 
nitrogen  and  20  percent  molecular  oxygen.  The  elemental  mass  fractions  are 

3.  5 (1-K„) 

= 4.5 


4.  5 


(46) 


K 


I 


0 


The  elemental  carbon  mass  fraction  is  the  parameter  chosen  to  charac- 
terize the  cases  and  is  varied  over  the  range 


O < < 0.  50 


(47) 


Allison,  D.O.,  "Calculation  of  Thermodynamic  Properties  of  Arbitrary 
Gas  Mixture  with  Modified  Vibrational-Rotational  Corrections.  " NASA 
TND.3538,  1966. 

Turner,  E.,  "The  NEST  Chemistry  Computer  Program,  " The  Aerospace 
Corporation,  Report  No.  TR-0059(6240-20)- 1,  Vol.  I,  July  1970. 
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This  system  will  range  from  oxygen  rich  (for  small  values  of  to 
carbon-rich  (for  large  K^).  The  dividing  line  between  oxygen-rich  and  car- 
bon-rich (with  respect  to  CO)  is  defined  from  Eq.  (2)  with  6 = 0.  This  yields 
K„  = 1/7  as  the  boundary  between  oxygen-rich  and  carbon-rich. 

The  first  set  of  comparisons  is  for  a carbon  mass  fraction  of  50  per-  | 

cent  and  is  contained  in  Figs,  4 through  6.  Figures  4 and  5 show  the  mixture  i 

molecular  weight  and  enthalpy  for  a range  of  temperature  and  pressure.  The 
error  at  low  pressure  and  high  temperatures  is  due  to  the  neglect  of  ioniza-  j 

tion.  The  vertical  slashes  on  the  curves  indicate  the  temperature  boundary  | 

below  which  condensed  carbon  will  form  (as  determined  from  the  exact  solu-  ! 

I 

tion  of  Ref.  8).  Thus,  the  approximate  solution  is  invalid  to  the  left  of  these 
slash  marks.  The  sudden  rise  in  enthalpy  in  this  regime  is  due  to  absorption 
of  the  excess  carbon  in  which  has  substantial  internal  energy  in  the  rota- 
tional and  vibrational  modes.  The  largest  error  is  about  three  percent,  and 
generally  the  error  is  less  than  one  percent. 

The  species  mass  fractions  at  100  atmospheres  are  shown  in  Figs.  6a 
and  6b.  In  the  range  where  solid  carbon  forms,  the  values  for  CO  and  C^  are 
obviously  wrong  (they  should  vanish  with  the  temperature).  Although  the 
error  in  the  species  mass  fractions  is  much  greater  than  for  the  mixture 
properties,  the  agreement  is  still  quite  good.  For  the  dominant  species 
(mass  fraction  greater  than  one-tenth)  the  largest  error  is  about  20  percent. 

The  next  set  of  comparisons  shows  the  variation  in  mixture  properties 
and  composition  with  carbon  mass  fraction.  Figures  7 and  8 show  the 
molecular  weight  and  enthalpy  for  two  pressures  and  several  temperatures. 

For  the  indicated  conditions,  neither  ionized  nor  condensed  species  are 
present,  and  the  agreement  is  excellent.  At  high  pressure  and  intermediate 
temperature,  slight  errors  are  evident  in  the  molecular  weight.  For  small 
values  of  carbon  mass  fraction,  the  error  is  due  to  the  neglect  of  nitric  oxide. 

The  error  at  large  carbon  mass  fractions  is  due  to  the  neglect  of  the  heavy 
molecular  weight  carbon-nitrogen  species  C^N^. 


P = 100  atm 


O EXACT  SOLUTION 
REF.  8 


r T = 10,000°K 


The  mixture  composition  is  presented  in  Figs.  9a  through  9d  for  a 
temperature  of  5000°K  and  pressures  of  1.0  and  100.0  atmospheres.  Again, 
the  largest  error  for  the  dominant  species  is  about  20  percent.  The  error 
may  be  much  larger  for  small  species  mass  fraction  (less  than  one-tenth), 
but  these  exert  no  appreciable  influence  on  the  mixture  properties. 

It  is  apparent  that  the  errors  are  smaller  for  the  lower  pressure.  This 
is  because  the  assumptions  and  approximations  utilized  in  the  model  are 
inherently  more  accurate  at  low  pressures.  That  is,  the  neglect  of  nitric 
oxide  and  is  an  increasingly  good  approximation  as  the  pressure  de- 

creases. The  ratio  of  cyano  to  cyanogen  will  increase  as  the  pressure 
decreases.  This  improves  the  accuracy  because  the  calculation  for  cyano  is 
more  accurate  than  that  for  cyanogen.  Finally,  the  presumption  that  the 
reactions  may  be  uncoupled  from  one  another  and  treated  separately  becomes 
a progressively  better  approximation  as  the  pressure  is  lowered. 

\ Figures  9a  through  9d  demonstrate  quite  clearly  that  the  partitioning  of 

1 the  problem  into  two  distinct  physical  regimes  (denoted  as  oxygen-  or  carbon- 

! rich)  is  an  excellent  approximation.  It  is  obvious  that  the  character  of  the 

I solution  changes  completely  as  approaches  1/7  (the  boundary  between 

I oxygen-  and  carbon-rich.  As  the  boundary  is  approached  from  small  values 

'■  of  K„,  the  oxygen  is  absorbed  into  CO,  and  the  free  oxygen  vanishes.  As 

■ the  boundary  is  approached  from  large  values  of  K^,  the  carbon-nitrogen 

j and  carbon  sublimation  species  vanish,  and  the  carbon  is  bound  up  in  CO.  It 

was  also  assumed  that  the  carbon  sublimation  species  would  vanish  faster 
^ than  the  carbon-nitrogen  compounds  as  the  oxygen- rich  boundary  was  ap- 

' proached.  This  is  indeed  seen  to  be  the  case.  Finally,  it  was  assumed  that 

l 

^ CO  would  always  form  to  the  maximum  possible  extent,  subject  to  the  availa- 

bility of  carbon  and  oxygen.  At  the  boundary  between  the  oxygen-  and  carbon- 
rich  domains,  the  formation  of  CO  is  unconstrained  by  the  available  carbon  or 
’ oxygen  and  will  reach  its  maximum  value  (for  a given  pressure  and  tempera- 

ture). It  is  apparent  from  Figs.  9a  and  9c  that  CO  does  indeed  reach  its 
maximum  value  at  € = 0. 
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The  possibility  of  ignoring  the  carbon-nitrogen  reactions  has  been 
alluded  to  several  times  in  Sections  II  and  III.  The  rationale  for  ignoring 
these  reactions  is  contained  in  Figs.  7 and  8.  The  solid  lines  include 
the  carbon-nitrogen  reactions  and  the  dashed  lines  ignore  them.  Where  only 
solid  lines  are  indicated,  the  difference  between  the  two  calculations  is  too 
small  to  distinguish  on  the  plot  scale  of  the  figures.  It  is  obvious  that  the 
differences  are  insignificant,  although  inclusion  of  the  reactions  does  slightly 
improve  the  accuracy. 

The  most  sensitive  test  of  ignoring  the  carbon-nitrogen  reactions  is 
encountered  for  a binary  system,  consisting  of  these  two  elements,  at  high 
pressure  and  intermediate  temperatures  (P  = 100  atm  and  T = 5000-6000°K). 

Figure  10  shows  such  a comparison  for  the  molecular  weight.  The  enthalpy 
shows  similar  errors.  Therefore,  if  only  the  mixture  properties  are  of 
interest,  the  carbon-nitrogen  reactions  may  be  ignored. 

B.  COMPUTATION  SPEED 

The  prime  motivation  of  this  report  was  to  develop  a computationally 
fast  (relative  to  exact  equilibrium  chemistry  codes)  method  for  the  mixture 
properties  of  a carbon-nitrogen-oxygen  system.  The  degree  to  which  this 
objective  was  met  is  examined  in  this  subsection. 

Particular  care  was  taken  to  make  the  comparisons  with  the  exact 
solutions  on  a one-to-one  basis.  Since  the  approximate  method  does  not 
consider  either  condensed  or  ionized  species,  these  were  also  excluded  from 
the  exact  solutions.  The  exact  solution  does  include  a large  number  of 
species  (20)  which  were  not  considered  in  the  approximate  model,  since  it  is 
precisely  the  inclusion  of  these  species  which  distinguishes  the  exact  solution 
from  the  approximate  solution.  Both  the  approximate  and  exact  solutions 
were  run  on  the  same  computer  (CDC  7600)  and  used  the  same  compiler 
(FTN).  In  both  cases,  only  the  actual  execution  time  was  included  (i.  e., 
input  and  output  time  was  excluded). 
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There  is  one  major  difference  between  the  exact  and  approximate 
solutions  which  has  a first-order  effect  on  the  time  comparisons  and  which 
cannot  be  treated  equivalently  in  both  calculations.  This  is  the  influence  of 
the  initial  guess  upon  the  computation  time  for  the  exact  solution.  Since  the 
approximate  solution  does  not  require  an  initial  guess,  there  is  no  way  to 
include  its  effect  upon  the  computation  time.  That  is,  even  if  the  approximate 
method  is  supplied  with  a good  initial  guess  for  the  solution,  there  is  no  pro- 
vision in  the  model  to  utilize  this  information.  The  exact  solution,  on  the 
other  hand,  is  strongly  dependent  upon  an  initial  guess.  Indeed,  if  the  initial 
guess  is  not  reasonably  accurate,  the  exact  solution  may  fail  to  converge. 

It  follows  that  the  computation  time  for  the  exact  solution  will  decrease 
dramatically  as  the  initial  guess  becomes  progressively  more  accurate.  This 
will  be  demonstrated  later  (see  Table  1). 

The  computation  time  for  the  approximate  solution  does  not  depend  on 
pressure  or  temperature.  It  is  also  independent  of  the  elemental  mass 
fractions,  unless  the  boundary  £ = 0 is  crossed.  That  is,  if  the  case  is 
carbon-rich,  the  cubic  which  yields  the  and  CO^  mass  fractions  is  not 
solved,  whereas  it  is  solved  for  the  oxygen- rich  case.  The  computation  time 
also  clearly  depends  on  whether  the  carbon-nitrogen  reactions  are  included. 
Since  the  mixture  properties  were  shown  to  be  essentially  unaffected  by  these 
reactions,  they  are  not  included  in  the  computation  time  for  the  approximate 
model.  The  pressure  equilibrium  constants  and  the  enthalpy  of  each  species 
were  treated  as  input  and  not  included  in  the  computation  time.  With  these 
provisos,  the  computation  times  for  a single  case  (given  pressure,  tempera- 
ture and  elemental  mass  fractions)  for  the  approximate  model  are  given  as 
follows : 


€ > 0:  Time  = 0.00021  sec 
f < 0:  Time  = 0.00016  sec 


(48) 
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It  is  emphasized  that  these  computation  times  are  quite  conservative, 
since  they  were  obtained  from  the  development  computer  code  which  was  not 
optimized  for  execution  time.  If  the  code  is  cleaned  up,  and  if  only  the  roots 
of  interest  are  computed  and  the  FORTRAN  coding  is  optimized,  the  compu- 
tation time  quoted  above  could  probably  be  reduced  by  a factor  of  two. 

The  exact  calculations  were  obtained  from  the  NEST  (Ref.  8)  chemistry 
code.  The  computation  times  were  not  strongly  dependent  upon  pressure, 
temperature  or  elemental  mass  fractions,  and  a single  condition  was  chosen 
for  comparative  purposes.  The  comparison  is  presented  for  a carbon-air 
mixture  at  a pressure  of  100  atm  and  a temperature  of  5000°K,  with  a 50 
percent  carbon  elemental  mass  fraction. 

As  discussed  above,  the  major  factor  influencing  the  computation  time 
is  the  accuracy  of  the  initial  guess.  Therefore,  the  equilibrium  composi- 
tion for  the  condition  of  interest  was  calculated  for  various  initial  guesses. 
The  initial  guesses  were  exact  solutions  for  conditions  close  to  the  condition 
of  interest.  The  computation  times  for  the  indicated  case  are  presented  in 
Table  1 as  a function  of  the  initial  guess. 

The  first  entry  in  the  table  uses  a preprogrammed  initial  guess  built 
into  the  NEST  code.  The  next  three  entries  use  initial  guesses  which  are 
exact  solutions  for  temperatures  which  differ  from  5000°K  by  varying 
amounts.  These  initial  guesses  have  the  same  carbon  mass  fraction  as  the 
t case  of  interest.  As  indicated,  the  computation  time  is  reduced  significantly 

^ as  the  initial  guess  is  improved.  The  last  two  cases  use  initial  conditions 

; where  both  the  temperature  and  carbon  mass  fraction  are  perturbed  about  the 

values  for  the  case  of  interest.  It  is  seen  that  a 10  percent  change  in  carbon 
mass  fraction  has  a significant  effect  upon  the  computation  time.  The 
initial  guesses  in  the  table  exhibit  the  kind  of  variations  encountered  during 
iteration  at  a boundary  layer  mesh  point.  A comparison  of  the  results  of 
Eq.  (48)  with  Table  1 indicates  the  analytic  solution  of  this  report  will  reduce 
I ' the  computation  time  by  two  or  three  orders  of  magnitude,  depending  upon  the 

accuracy  of  the  initial  guess  available  to  the  exact  solution.  For  typical 
! chemically  reacting  boundary  layer  calculations,  a reduction  in  computation 

time  of  the  order  of  500  is  a reasonable  expectation. 


Table  1.  Computation  Time  for  Exact  Solution  of 
Carbon-Air  Gas  Mixture 

P =100  atm 

T = 5000°K 

= 0.50 


Initial  Guess  Gonditions 

Computation  Time  (sec) 
for  Indicated  Case 

Not  Applicable 

0.295 

T = 6000  = 0.50 

0.  032 

T = 5500  = 0.50 

0.027 

T = 5100  = 0.50 

0.  on 

T = 5500  = 0.45 

0.  157 

T = 5100  = 0.45 

0.043 

SECTION  V 


SUMMARY  AND  CONCLUSIONS 

An  approximate  analytic  solution  has  been  developed  for  the  equilibrium 
composition  of  a carbon-nitrogen- oxygen  gas  mixture  (including  provision  for 
an  inert  monatomic  species).  This  was  accomplished  by  partitioning  the 
problem  into  carbon- rich  and  oxygen- rich  domains  and  then  subdividing  each 
domain  into  high  and  low  temperature  regimes.  In  this  manner,  the  dominant 
reactions  are  uncoupled  from  one  another  and  treated  sequentially  in  pairs. 
This  allows  the  solution  to  be  expressed  in  terms  of  cubics  and  quartics, 
both,  of  which  possess  analytic  solutions. 

Carbon-nitrogen  reactions  were  included  in  the  formulation  and  were 
found  to  exert  a negligible  influence  on  the  mixture  properties.  Ignoring 
these  reactions  yields  the  solution  in  its  simplest  form,  namely,  two  cubics 
and  a quartic  (for  pressure  and  temperature  as  independent  variables).  The 
solutions  are  trivial,  and  there  is  no  ambiguity  concerning  which  is  the 
desired  root. 

The  formulation  was  carried  out  for  two  different  sets  of  independent 
variables:  pressure  and  temperature,  and  density  and  temperature.  The 
use  of  density  as  a variable  uncoupled  the  equilibrium  relationships  from  the 
molecular  weight  and  substantially  reduced  the  algebraic  complexity.  The 
uncoupling  allowed  the  solution  for  the  N2  and  CO  reactions  to  be  obtained  in 
terms  of  two  uncoupled  quadratics  instead  of  the  quartic  which  resulted  from 
the  use  of  pressure  as  an  independent  variable.  With  density  as  a variable, 
the  solution  (ignoring  carbon-nitrogen  reactions)  was  obtained  in  terms  of 
two  cubics  and  two  quadratics. 

Extensive  comparisons  were  made  with  exact  solutions  to  determine 
the  accuracy  of  the  approximate  solution.  The  maximum  error  in  mixture 
properties  (molecular  weight  and  enthalpy)  was  about  four  percent  (ignoring 
carbon-nitrogen  reactions),  and  generally  the  error  is  less  than  one  percent. 


L 
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Since  the  required  thermodynamic  data  (heats  of  formation,  vibrational  and 
rotational  constants,  etc.  ) are  not  known  to  within  five  or  10  percent,  the 
difference  between  the  exact  and  approximate  solutions  is  less  than  the  un- 
certainty in  the  input  data.  In  this  context,  the  approximate  solution  may  be 
regarded  as  exact. 

A detailed  examination  was  made  of  the  computation  times  for  the 
approximate  and  exact  solutions.  It  was  found  that  the  analytic  solution  re- 
duces the  computation  time  between  two  and  three  orders  of  magnitude, 
depending  upon  the  accuracy  of  the  initial  guess  available  to  the  exact  solu- 
tion. For  a boundary  layer  application,  the  analytic  solution  will  reduce  the 
computation  time  required  to  determine  the  mixture  properties  by  approxi- 
mately a factor  of  500. 

In  conclusion,  it  is  believed  that  the  analytic  solution  of  this  report 
represents  a significant  advance  in  the  ability  to  handle  chemically  reacting 
flows  (for  a carbon-nitrogen-oxygen-inert  system)  in  a computationally 
efficient  manner,  with  essentially  no  loss  in  accuracy.  For  fluid  dynamics 
problems  of  fixed  elemental  composition,  the  computation  speed  is  probably 
competitive  with  the  two-dimensional  interpolation  required  for  use  with 
precomputed  Mollier  diagrams.  For  this  application,  the  storage  require- 
ments for  the  Mollier  diagram  may  be  eliminated.  For  problems  of  varying 
elemental  composition,  the  reduction  in  calculation  time  is  so  great  (com- 
t pared  to  exact  solutions)  that  the  use  of  the  analytic  solution  essentially 

eliminates  computational  constraints  imposed  by  the  chemistry. 


L. 


APPENDIX  A 


PRESSURE  EQUILIBRIUM  CONSTANTS 


The  reactions  considered  are  all  of  the  following  form: 


(A-1) 


where  the  A^  are  the  reactants  and  the  products.  The  pressure  equilibrium 
constant  for  this  reaction  is  defined  in  terms  of  the  partial  pressures  as 


K 


P 


1 2 

P ^B^)  P “"(B^) 

^1  ^2 

P ^Aj)  P ^{A^) 


(A-2) 


The  partial  pressure  of  the  species  A^^  is  denoted  as  P(A^)  and  is  related  to 
the  species  mass  fraction  by  the  following: 


K. 

P(A.)  = MP  (A-3) 


where  P is  the  pressure  and  M is  the  mixture  molecular  weight. 

Substitution  of  Eq.  (A-3)  into  Eq.  (A-2)  yields  the  «.  coefficient  which 
is  defined  in  Section  II. 

The  equilibrium  constants  are  computed  from  the  partition  functions  as 
follows: 

-AE 

log  Kp  = + bjlog  Qp(Bj)  + b^  log  Qp(B2)  + 

- a^  log  Qp(Aj)  - a^  log  Qp(A2)  (A-4) 

where  Eq  denotes  the  heat  of  formation. 


The  partition  functions  are  given  in  Table  A-1  for  the  species  of 
interest.  The  logarithms  are  all  natural  logarithms  and  the  temperatures 
are  in  degrees  Kelvin.  The  thermodynamics  data  was  obtained  from  Refs.  7 
and  9.  For  the  diatomic  species,  the  partition  functions  were  simplified  by 
assuming  that  the  rotational  and  vibrational  partition  functions  for  each 
electronic  state  are  approximately  equal  to  the  ground  state  (Ref.  9).  Fur- 
thermore, all  the  expressions  were  truncated  at  a level  consistent  with  the 
temperature  range  of  interest  (0-10,000°K), 

The  values  of  AE^  are  tabulated  below  for  the  specific  reactions  of 
interest. 


Kp 

AE„ 

(°K) 

1 

67190 

2 

59350 

3 

128903 

4 

113228 

5 

72439 

6 

87305 

7 

72763 

8 

162005 

9 

75482 

9.  Horton,  T.E.,  "The  Computation  of  Partition  Functions  and  Thermo- 
chemistry Data  for  Atomic,  Ionic,  Diatomic  and  Polyatomic  Species, 
Jet  Propulsion  Laboratory,  TR  32-1425,  1970. 


Table  A-1.  Partition  Functions 


log  Qp(NO)  = 3.5  log  T + 0.55  + log[2  + 2 exp(-174/T)]  - log[l  - exp(-271 9/T)] 


n 


^^cecliA'i^  ^ar/T/f  ~ 


APPENDIX  B 


SOLUTION  FOR  CO^  AND  O^  REACTIONS 


The  solution  for  the  equation  set  of  Subsection  II-A-I  is  presented  in 
this  appendix.  The  following  variables  are  defined: 


X = MK 


CO 


(B-I) 


Y s mK 


O 


Substitution  of  these  variables  into  Eqs.  (4)  through  (6)  yields  a cubic  in  Y 

AY^  + BY^  + CY  + -^6  = 0 


C 8 O 

a, 


B.-TT^(Kc-X-£)-Tf 


(B-2) 


After  the  cubic  is  solved,  X is  obtained  from 


X = 


7'  t-iT’f'^^(Ko-|Kc) 


(B-3) 
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The  mass  fractions  for  CO  and  O are  calculated  from  Eq.  (B-1)  and 
the  mass  fractions  of  and  CO^  from  Eq.  (6)  in  Subsection  II-A-1. 


The  asymptotic  behavior  of  the  above  system,  for  large  and  small 
temperatures,  is  obtained  from  a knowledge  of  the  behavior  of  the  01  coeffi- 
cients. That  is,  as  the  temperature  vanishes,  the  Ct  coefficients  become 
infinite,  and  as  the  temperature  becomes  infinite,  the  01  coefficients  vanish. 

Therefore,  as  the  temperature  becomes  infinite,  Eq.  (B-2)  yields 

As  the  temperature  vanishes,  Y vanishes  and  the  cubic  and  linear 
terms  of  Eq.  (B-2)  become  predominant.  Therefore,  the  cubic  may  be  ) 


solved  to  yield 


T-*-0:  Ot^Y^ 


-l^o) 


(B-6) 


i Substitution  of  Eqs.  (B-6)  into  (B-3)  shows  that  X also  vanishes  as  the 

1 temperature  does.  Then,  using  Eqs.  (B-5)  and  (B-6),  the  asymptotic 

f behavior  of  Eqs.  (7)  and  (8)  is  obtained. 

I 

: As  £ approaches  zero,  it  is  obvious  that  Y must  also  vanish  and,  thus, 

[ 

the  behavior  of  Eq.  (9)  is  obtained. 


Some  comments  about  the  roots  of  the  cubic  are  in  order.  The  correct 
root  must  lie  in  the  range 


For  all  the  cases  considered  in  this  report,  there  was  only  one  root  in 
the  allowable  range.  Generally,  there  were  two  negative  roots  and  one  posi- 
tive root.  If,  for  some  case,  more  than  one  root  lies  in  the  allowable  range, 
the  correct  root  should  yield  a positive  value  for  X, 


[ 

i 


/v^r 


APPENDIX  C 


SOLUTION  FOR  N2  AND  CO  REACTIONS 


The  solution  for  the  set  of  Eqs.  (11)  in  Subsection  II-A-2  is  presented  in 
this  appendix.  The  following  variables  are  defined 


X = KqM 


Y s K^M 


Z 


The  equations  may  be  rearranged  to  yield  a quartic  in  Z. 


AZ^  + BZ^  + CZ^  + DZ  + E = 0 


A = a3a4  (n^  + + ^3) 

B = ct^a^[{rj^  + + rj^)  + (t)^  + + rj^)] 

C s a^{r)^  + 7)3) (7)j  + n^)  . + 773  + ^4)  + 

- a^a^iZT]^  + 773  + 77^) 

D =-a3(277^  + 773  + 77^)  . ^ (rj^  + + 77^) 


(C-1) 


(C-2) 


The  7}  coefficients  are  defined  by 


^ - 1 . I 
14  ~ 

KnMi 


1 I 


(C-3) 


After  solving  the  quartic  for  Z,  X and  Y are  calculated  from 


Y = 12  [1  - Z{r]^  + T]^)  - ct^'zriVz  + ^3)] 
X = 16  [1  - Z(77j  + n^)  - «4Z^(772  + 774)] 


(C-4) 


The  molecular  weight  is  calculated  from 


= 5^- (z 


(C-5) 


Equations  (C-1)  are  used  to  calculate  the  mass  fractions  for  atomic 
oxygen,  carbon  and  nitrogen.  The  remaining  mass  fractions  are  calculated 
from  Eq.  (11). 
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Noting  that  the  Ot  coefficients  are  infinite  for  zero  temperature  and  zero 
for  infinite  temperature,  the  asymptotic  behavior  of  the  equations  is  examined. 
As  the  temperature  vanishes,  the  even  power  terms  of  the  quartic  are  dominant 
and  the  equation  may  be  recast  as  a quadratic.  That  is,  define  a variable  F as 

F=  (C-6) 

Substitution  of  this  expression  into  the  quartic  yields,  for  vanishing  tempera- 
ture 

^3^  “ ^<^^2  + + r]^)  + I = 0 (C-7) 

The  roots  are 


F = 


’ Vo  + V. 


(C.8) 


The  first  root  corresponds  to  the  oxygen-rich  case  and  the  second  to 
the  carbon-rich  case.  Substitution  of  Eq.  (C-8)  into  Eq.  (C-4)  yields 

t >0  (oxygen- rich) : 


Z = 0 

“4^'  = 
Y = 0 
X = 16 


V2  +^3 


^3  - ^4 
^2  + ^3 


.73- 


(C-9  continued) 


€ < 0 (carbon-rich); 


Z = 0 


a.z^  = 1 


'4  r)^  + 774 

X = 0 

774  - ^3 


Y = 12 


\ + % 


The  results  of  Eq.  (12)  are  easily  obtained  from  Eq.  (C-9)  for  vanishing 
temperature. 

As  the  temperature  approaches  infinity,  the  quartic  yields 

1 


Z = 


77i  + 773  + 774 

This  result  is  used  in  Eqs.  (C-4)  and  (C-5)  to  obtain  the  value  of  Eq.  (13), 
The  correct  root  of  the  quartic  lies  in  the  range 


(C-10) 


0 < z < 


>71  +773  + 774 


(C-11) 


Since  an  analytic  solution  for  a quartic  exists,  there  is  no  difficulty  in 
obtaining  the  four  roots.  The  analytic  solution  is  only  slightly  more  complex 
than  that  of  a cubic  and  requires  negligible  computation  time  to  obtain  the 
roots.  At  high  temperature,  only  one  root  lies  in  the  allowable  range.  How- 
ever, at  low  to  moderate  temperatures,  there  are  two  roots  in  the  allowable 
range  (one  for  €>0  and  one  for  € <0).  The  correct  root  is  that  value  which 
yields  positive  values  for  X and  Y.  For  the  cases  considered  in  this  report, 
there  were  generally  two  positive  and  two  negative  roots.  When  two  positive 
roots  are  in  the  allowable  range  of  Eq.  (C-11),  it  is  obvious  from  Eqs.  (C-4) 
that  the  correct  root  is  the  smaller  of  the  two  values. 
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0: 
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Solution  for  K., 

N 

It  is  apparent  that  the  above  equal  ions  are  singular  if  the  nitrogen 
elemental  mass  fraction  is  zero.  Thus,  for  = 0,  the  CO  dissociation 
reaction  is  considered  separately.  The  solution  is  obtained  in  terms  of  the 
following  quadratic: 


(C-12) 


At  low  temperatures,  the  plus  sign  is  used  with  the  quadratic  solution  radical 
for  the  oxygen-rich  branch  and  the  minus  sign  for  the  carbon- rich  solution. 

At  high  temperatures  only  one  positive  root  exists,  and  it  is  the  solution. 

The  molecular  weight  is  calculated  from 

M (C-13) 

Mj 


I • 


i 


Y.K^M 


(D-1) 


Z s 


The  solution  is  obtained  in  terms  of  a quartic 
AZ'^  + BZ^  + CZ^  + DZ  + E = 0 


A = ^ a^alrjl 

B = 2 ^ - 13^4^6 ''■  ^1^ 

_ 2 “4 

C = - 2 a4a5r74  - (^3  + ^74)  + 


(D-2) 


J!  jjWfcfejM!  y.  ■■■ 


The  r]  coefficients  are  defined  by 


1 h 

77  s j.  - ^ - x ZL 

A * ^ X £ TT- 


KnMi  Kj 


\ -n 


24  K., 


{D.3) 


^74=  ^2  + 


775-  /7i+^3 


After  the  quartic  is  solved  for  Z,  the  following  expressions  are  used  to 


calculate  Y and  M. 


„ 24(i-,5Z.,4«^z2) 

, . 168  „ ^ „ 

1 13  ^4^6^ 

M = g^-fz+a4z2+^a(,Yzl 

K«-  L J 


(D.4) 


After  Eqs.  (D-1)  are  used  to  obtain  the  mass  fractions  of  atomic  nitrogen 
and  carbon,  the  remaining  mass  fractions  are  obtained  from  Eqs.  (15). 


The  asymptotic  behavior  of  the  system,  at  low  temperatures,  is  again 
obtained  from  the  knowledge  that  the  a coefficients  become  infinite  as  the 
temperature  vanishes.  It  is  also  necessary  to  know  the  relative  rates  at 
which  the  Cl  coefficients  become  infinite.  The  results  of  Appendix  A are  used 
to  obtain  the  following  expression; 


►0:  — — = 0 


(D-5) 


When  the  substitution  F = ^4^  used,  the  quartic  becomes  a quadratic 
at  low  temperatures  and  may  be  solved  directly  for  F.  Both  roots  of  the 
quadratic  are  equal  and  yield 

F = (D-6) 

''4 

Thus,  in  contrast  to  the  results  of  Appendix  C for  the  N2  and  CO  reactions, 
only  one  valid  solution  branch  exists  at  low  temperatures. 

Using  the  result  for  F yields,  as  the  temperature  vanishes 


Y = 0 


oe^YZ  = 0 


(D-7) 


As  the  temperature  becomes  infinite,  the  quartic  is  linear  in  Z and 


y ields 


^3  +^5 


(D.8) 


With  the  above  expressions,  the  limiting  behavior  of  Eqs.  (17)  in 
Subsection  II-B-1  is  easily  obtained. 

As  before,  the  correct  root  of  the  quartic  lies  between  the  values  for 
the  two  temperature  extremes.  That  is 


O 


<Z  < 


1 

^3+^5 


(D-9) 


For  the  cases  considered  in  this  report,  there  were  generally  two  posi- 
tive roots  in  the  allowable  range.  Again,  the  correct  root  yields  a positive 
value  for  Y.  From  Eq.  {D-4)  it  is  apparent  that  thi  correct  root  will  be  the 
smaller  of  the  two  roots  in  the  allowable  range. 
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APPENDIX  E 


SOLUTION  FOR  AND  REACTIONS 


The  solution  of  Eqs.  (22)  in  Subsection  II-B-3  is  obtained  in  this 
apjjendix.  The  transformed  variable  is  defined  by 


Y H K^M 


(E-1) 


The  equations  yield  the  following  cubic: 


AY^  + BY^  + CY  + D = 0 


/ K K - 

A = OCg  ^^^5  + j;;q  + 35 


^I  - ^4 


(E-2) 


Kf 

^ ^5  Mj  ^ 12 


D. 


The  molecular  weight  is  obtained  from 


M = 


^C-^4 


Y +a^Y^  +agY^ 


(E-3) 


After  solving  the  above  equations  for  the  atomic  carbon  mass  fraction, 
the  mass  fractions  for  and  are  calculated  from  the  equilibrium  rela- 
tionships of  Eqs.  (22). 
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As  the  temperature  becomes  infinite,  the  cubic  yields 

Y = ^ (E-7) 

With  these  results,  the  behavior  of  Eq.  (24)  is  easily  obtained. 

The  only  rem.aining  limit  of  .nterest  occurs  for  the  situation  where  0. 

approaches  K„.  This  indicates  t]iat  all  the  carbon  is  contained  in  CO  and  the 

carbon-nitrogen  compounds  and  is  not  available  for  the  formation  of  and 

C-.  Thus,  foi  8.  ^ K„,  the  mass  fractions  for  and  C-  are  zeroed  and 
3 ' ^^4  C 2 3 

the  cubic  is  not  solved. 
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The  correct  root  of  the  cubic  lies  in  the  range  defined  by 


0£Y<^ 


(E-8) 


At  low  temperatures  there  is  only  one  real  root,  whereas  at  high 
temperatures  there  are  three  real  roots,  of  which  only  one  lies  in  the  allow- 
able range. 
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APPENDIX  F 

NITROGEN-OXYGEN  SYSTEM 


1 
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As  mentioned  earlier,  the  analysis  of  Sections  II  and  III  yields  an 
approximate  model  for  the  nitrogen- oxygen  system  if  the  carbon  elemental 
mass  fraction  is  zero.  If  pressure  is  the  independent  variable,  the  approxi- 
mations consist  of  ignoring  the  formation  of  nitric  oxide  and  uncoupling  the 
N^  and  O^  dissociation  reactions.  If  density  is  the  variable,  no  approxima- 
tions are  required  to  uncouple  the  reactions.  Again,  nitric  oxide  is  ignored. 

These  approximations  have  been  used  by  several  investigators  to 
obtain  analytic  solutions  for  the  mixture  properties  of  air  with  generally 
good  results.  There  are,  however,  situations  in  which  these  assumptions 
begin  to  influence  the  results  (particularly  at  high  pressures).  It  is  the 
purpose  of  this  appendix  to  determine  the  equilibrium  composition  of  a 
nitrogen-oxygen  system  in  as  general  a manner  as  possible  within  the  con- 
straint of  obtaining  an  analytic  solution.  As  before,  this  restriction  requires 
the  solution  to  be  expressed  in  terms  of  a polynomial  of  fourth  order  or  less. 
The  formulation  is  first  carried  out  using  pressure  and  temperature  as  the 
independent  variables.  Following  this,  the  problem  is  reformulated  using 
density  and  temperature  as  the  independent  variables. 

1.  PRESSURE  AND  TEMPERATURE  AS  INDEPENDENT 

VARIABLES 

Hansen  (Ref,  1)  and  Mikhailov  (Ref.  3)  ignored  the  formation  of  nitric 
oxide  and  treated  the  oxygen  and  nitrogen  dissociation  reactions  sequentially 
to  obtain  the  properties  of  air.  No  provision  was  made  for  the  inclusion  of  an 
inert  species.  The  solution  was  obtained  in  the  form  of  quadratic  equations. 
Although  Hansen's  and  Mikhailov's  results  are  specifically  for  air,  they  are 
easily  generalized  to  other  nitrogen-oxygen  ratios. 


i; 
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The  present  solution  includes  nitric  oxide,  treats  all  the  reactions 
simultaneously  and  includes  an  inert  species.  To  the  author's  knowledge, 
the  present  solution  is  the  most  general  analytic  solution  available  for  a 
nitrogen-oxygen  system.  The  penalty  for  this  generality  is  that  the  solution 
is  obtained  in  terms  of  a quartic  equation  instead  of  a quadratic.  Since 
Hansen's  and  Mikhailov's  results  are  quite  accurate  for  most  applications, 
the  present  solution  would  be  used  only  at  high  pressures  (greater  than  100 
atmospheres)  or  if  one  were  interested  in  the  NO  mass  fraction. 


The  elemental  mass  balances  are  given  by 


o 

= 

+ 

+ 1^ 

K ^ 

O 

^2 

30 

NO 

= K 

+ 

K ^ 

N 

N 

^2 

30 

NO 

(F-1) 


Ki  =Kj 


The  elemental  mass  fractions  (denoted  by  a superscript  tilde)  are 
presumed  known  and  must  sum  to  unity. 


The  equilibrium  relationships  are 


= o^mkI 

"^2 

2 O 

= 

"^2 

4 N 

= a„M 


(F-2) 
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The  other  Cl  coefficients  are  defined  in  the  main  text.  The  pressure 


equilibrium  constant  Kp  is  defined  for  the  reaction  given  by 

9 


NO-^N  + O 


(F-3) 


The  molecular  weight  is  given  by 
K 


J_  “O  

M ' 32 


''°2  . . ■'no  , ■'l 


14  28 


30  M 


(F-4) 
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Transformed  variables  are  defined  by 


X = MK 


O 


(F-S) 


Y s MK 


N 


The  above  equations  may  be  combined  to  yield  the  solution  in  the  form 
of  a quartic. 


AX^  + BX^  + CX^  + DX  + E = 0 


A 2 56  2 

A = ci^a^  - doOr 


'4''2  225  “2^9 


B _ 2 225  ^^2  ' + 


15“2“9^^  + Xj  + X^) 


C s 


8 

15 


«9  X2(Xi  - 1)  - (Xj  + D^j 


225  “9  ^^1  “ ^2^  " 7 “2^1^^  ^2^ 


- 64012014X1  +«4(1  + Xj)' 


(F-6  continued) 
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A . ^ ■ J . v" 


D = -^  Xjttgd  + Xj  - ZX^)  - 6404X^(1  +Xj)  + 

- -y  Xi  (1  + Xj  + X^) 

E = 256  Xj  (4a^  + y) 


The  X coefficients  are  given  by 

#«v 

K, 


2Ki 


'O 


(F-7) 


K 


N 


14;5 


After  solving  the  quartic  for  X,  the  transformed  variable  Y and  molecu- 
lar weight  are  obtained  from 


(F-8) 


+y(i  tX2+^a,x)  (^-2)  - 0 

The  limiting  behavior,  for  large  and  small  temperatures,  is  given  by 


2 2K 

T_0:  * -f- 


32  X, 


(F-9) 


T-*.oo:  X 


I + Xi  + X2 
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From  these  results,  it  is  easily  seen  that  the  mass  fraction  of  nitric 
oxide  vanishes  at  high  and  low  temperatures.  At  low  temperatures  the  nitro- 
gen and  oxygen  are  bound  up  in  molecular  form  and  are  contained  in  atomic 
form  at  high  temperatures. 

The  solution  for  X is  not  monotonic  between  the  temperature  extremes. 
Therefore,  an  upper  bound  for  X may  be  defined  by  using  the  maximum 
values  for  and  the  molecular  weight.  The  correct  root  lies  in  the  follow- 
ing range: 


O <X  < 2 

P 


(F-10) 


If  more  than  one  positive  root  lies  in  the  allowable  range,  the  correct 
value  will  be  the  smallest  root  (i.e.,  that  root  which  yields  a positive  value 
for  Y). 

2.  DENSITY  AND  TEMPERATURE  AS  INDEPENDENT 

VARIABLE^ 

Hochstim  (Ref.  2)  has  obtained  the  most  general  solution  for  this  set  of 
independent  variables.  The  solution  included  nitric  oxide  and  treated  the 
reactions  simultaneously.  Thus,  Hochstim's  solution  is  an  exact  analog  of 
that  obtained  above  for  given  pressure  and  temperature.  Since  the  molecular 
weight  is  not  required  for  the  determination  of  the  mass  fractions  (for 
specified  density  and  temperature),  the  inert  species  has  no  influence  on  the 
solution  until  the  molecular  weight  is  calculated. 

In  the  following,  the  only  difference  between  this  analysis  and  Hochstim's 
is  that  here  the  solution  is  specified  explicitly  in  terms  of  a quartic,  whereas 
his  solution  was  presented  as  two  coupled  quadratics  which  were  solved 
iteratively. 
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The  molecular  weight  in  the  equilibrium  relationships  of  Eq.  (F-2)  is 
absorbed  in  a new  set  of 3^^  given  by 


(F-11) 


“2  = ^"2 


ar^  = Ma4 


The  a contain  the  product  of  pressure  and  molecular  weight.  The  ideal 
gas  law  equates  this  product  to  the  product  of  density,  temperature  and 
universal  gas  constant.  Since  density  and  temperature  are  given,  the  product 
of  pressure  and  molecular  weight  is  known;  therefore,  the  3.  are  known 
quantities . 

Use  of  the  3^  in  the  equilibrium  relationships  and  substitution  into 
Eqs.  (F-1)  yields  the  solution  in  the  form  of  a quartic. 

AK^  + BK^  + CK^  + DKq  + = 0 

B s 2 - 225^9  ~ '^^“9 


C3  a4-20^a4Ko+^a2Ko  + 

-L*  64  g2  ~ 

■ 15  “9  ■ 225^9  N 


D.  ^a^K^.2a4K^ 


(F-12) 


-90- 


r- 


After  the  quartic  for  is  solved,  the  nitrogen  mass  fraction  is 


obtained  from 


4 N N 


r 7 ^ 


0 


(F-13) 


RESULTS 


The  analytic  solution  of  Subsection  F-1  was  compared  with  an  exact 
solution  (Ref.  8)  for  argon  free  air.  The  air  was  composed  of  20  percent  O, 
and  80  percent  N2  by  volume.  This  corresponds  to  the  following  mass 
fractions : 


Kq  = 0.2222 


= 0.7778 


The  species  mass  fractions  are  compared  in  Fig.  F-1,  where  the  solid 
lines  are  the  solution  of  Subsection  F-1  and  the  symbols  denote  the  exact 
solution.  The  agreement  is  seen  to  be  excellent. 
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